Multifractal dimensions for all moments for certain critical random matrix ensembles 

in the strong multifractality regime 
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We construct perturbation series for the qth moment of eigenfunctions of various critical random 
matrix ensembles in the strong multifractality regime close to localization. Contrary to previous in- 
vestigations, our results are valid in the region q < 1/2. Our findings allow to verify, at first leading 
orders in the strong multifractality limit, the symmetry relation for anomalous fractal dimensions 
A g = Ai_ ? , recently conjectured for critical models where an analogue of the metal-insulator transi- 
tion takes place. It is known that this relation is verified at leading order in the weak multifractality 
regime. Our results thus indicate that this symmetry holds in both limits of small and large coupling 
constant. For general values of the coupling constant we present careful numerical verifications of 
this symmetry relation for different critical random matrix ensembles. We also present an example 
of a system closely related to one of these critical ensembles, but where the symmetry relation, at 
least numerically, is not fulfilled. 

PACS numbers: 71.30.+h, 05.45.-a, 05.45.Df 



I. INTRODUCTION 

It was realized more than fifty years ago by Anderson 
[l| that disorder can lead to localization of wave func- 
tions in quantum systems. Depending on the strength 
of disorder, three-dimensional wave functions at a given 
energy are either all delocalized (metallic phase) or all 
localized (insulator phase). In recent years, such An- 
derson transitions have been observed experimentally in 
various media: localization of Bose-Einstein condensates 
into disorder created by laser speckle [2j, localization of 
ultrasound in elastic media [|| , localization of cold atoms 
in a kicked potential 

Many works have investigated theoretically the proper- 
ties of quantum disordered systems exactly at the metal- 
insulator transition A characteristic feature of wave 
functions at the transition point is that they are nei- 
ther localized nor delocalized. Their localization prop- 
erties are described by a set of multifractal exponents 
characterizing the behavior of the wave function mo- 
ments at different scales. More precisely, if we consider a 
wave function |W) belonging to an iV-dimensional Hilbert 
space with components 9j, then the qth moment of the 
wave function behaves for large N as Ylj=i \^j\ 2q rx 
N -D q (q-i)_ If |^ is i oca i ize H then Dq = for all pos- 
itive q, while if the state is delocalized over the whole 
Hilbert space then D q is the dimension of the system it- 
self. Multifractality is characterized by the fact that D q 
is a a nontrivial function of q. 

Understanding the properties of multifractal dimen- 
sions has been the subject of intense research in the past 
few years (see e.g. Q and references therein). In order 
to investigate the mechanisms involved, various models 
have been proposed. The simplest models where an ana- 
logue of the metal-insulator transition can be observed 
and studied are random Hermitian or unitary matrices 
whose off-diagonal elements M mn decrease as a power of 



the distance from the diagonal, as 
9 



m- n > 1, 



(1) 



where g plays the role of a coupling constant (here m, 
n are d-dimensional vectors, where d is the space di- 
mension). When r > d, eigenvectors of these matrices 
are localized; they are delocalized for r < d. Anderson- 
like transition occurs in the localization properties of the 
eigenvectors at r = d, and typically a multifractal behav- 
ior is observed Q, similar to the Anderson model at the 
point of metal- insulator transition 0]. Such ensembles 
for r — d are often called critical random matrix ensem- 
bles. The main attention in this field was given to the 
investigation of the ensemble of so-called critical banded 
random matrices (CrBRM) 8], whose matrix elements 
M mn are independently distributed Gaussian variables 
with zero mean and variance decreasing away from the 
main diagonal as in (|T|) with r = d = 1. Other mod- 
els with similar properties, such as ultrametric critical 
random matrix ensembles, have also been proposed 
Recently other critical random matrix ensembles, based 
on Lax matrices of classical one-dimensional iV-particle 
integrable systems, were introduced in (lpj . It was shown 
that the spectrum of such matrices is of intermediate 
type [ll|, and that the corresponding eigenvectors dis- 
play multifractal properties |T2| . 

Most of the analytical results for these critical ensem- 
bles are obtained within perturbation series expansions 
in the regime of strong multifractality where eigenvectors 
are almost localized, or in the regime of weak multifrac- 
tality where eigenvectors are almost entirely extended. 
Nevertheless, a few relations [HI, [l4| are believed to be 
universal and valid for a large (but not well defined) class 
of critical ensembles. One of the most striking ones in- 
volves anomalous fractal dimensions A q defined by 



A q = (D g -d)(q-l) 



(2) 
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which characterize the departure of multifractal expo- 
nents from their value in the delocalized case (since 
A q = for extended states). In pj| it was conjectured 
that these anomalous dimensions have to be symmetric 
with respect to the axis q = 1/2, namely 



Al-g. 



(3) 



One of the analytical arguments in favor of this result is 
the four-loop calculation of A q within nonlinear er-modcl 
in 2 + e dimensions [l5l ] , which manifests symmetry ([3]) . 
Recently, it has been argued that this relation holds ex- 
actly in models which can be mapped onto nonlinear a- 
models [l6j . This is for instance the case for critical ran- 
dom matrix models in the regime of weak multifractality. 
The fact that this symmetry is not always verified is il- 
lustrated by the existence of critical systems for which 
such a symmetry relation does not hold in general [17| 
(see also Section IVII|) . Thus the domain of validity of ([3]) 
remains unclear. 

For general systems, the perturbation expansion ap- 
proach remains the main analytical tool for investigating 
multifractal dimensions. For systems of the form (JT|), an- 
alytical expressions for multifractal dimensions at leading 
orders of perturbation expansion have been obtained in 
the domain q > 1 /2, both in the strong and weak multi- 
fractality regimes |12| . While the methods used to derive 
these results are still applicable in the domain q < 1/2 in 
the weak multifractality regime, no perturbation expan- 
sion approach has been developed (to the best of authors' 
knowledge) in the domain q < 1/2 in the strong multi- 
fractality case. Therefore direct analytical verification of 
the symmetry ([3]) in the regime of small g was not yet 
possible. 

The goal of this paper is to derive analytical expres- 
sions for the D q in the domain q < 1/2 where previous 
approaches did break down. We obtain a perturbative 
expression for q < 1/2 in the regime of strong multi- 
fractality, and perform extensive numerical calculations 
which corroborate our analytical findings. This enables 
us to test analytically the correctness of the symmetry 
relation ([3j for several different critical random matrix 
ensembles. The main conclusion of the paper is that, 
at least at leading orders in the perturbation series, the 
symmetry is fulfilled. 

The outline of the paper is the following. In section 
HIl the models investigated are introduced. Then mul- 
tifractal dimensions D q are defined in section IIIII The 
perturbation expansion of D q in the strong multifractal- 
ity regime for q < 1/2 yields a nontrivial result even at 
lowest order; the corresponding calculation is done in sec- 
tion [IVJ In section [V] the calculation is performed at the 
next-to-leading order. This enables to demonstrate the 
validity of the symmetry relation (j3|) up to this order for 
strong multifractality. Together with known results for 
the weak multifractality regime, this proves that the sym- 
metry holds in both extreme regimes. In order to assess 
the validity of this relation in general, numerical results 
are presented for the different models considered in sec- 



tion I VII In section IVIII we present numerical results for 
a system without perturbation theory expansion, where 
eigenfunctions display nontrivial multifractal dimensions, 
but for which the symmetry property ([3]) is not valid. In 
section [Villi a toy model related to the above ensembles, 
and where exact analytical calculation of multifractal di- 
mensions is possible, is briefly discussed. 



II. THE MODELS 

We now introduce the models that we consider in this 
paper. Four critical random matrix ensembles associ- 
ated with the Lax matrices of classical integrable systems 
have been introduced in [l2j. Calogero-Moser (CM) N- 
particle systems yield three ensembles of N x N Hermi- 
tian matrices of the form 



M„ 



■ PmSmn + - S rnn )V(m - 77,), 



(4) 



where p m are independent Gaussian random variables 
with zero mean and unit variance, and V(k) is one of the 
three following functions 



/' 



Nsinh(fik/N) 



Nsin(fik/N) 



(5) 



with [i a parameter (S mn is the Kronecker symbol). We 
denote these ensembles respectively by CM r , CM/j and 
CM t . The fourth ensemble introduced in [lj| corresponds 
to the Lax matrix of Ruijsenaars-Schneider (RS) system. 
It consists of a set of unitary N x N matrices which, up 
to a simple transformation, reduce to 



M„ 



sin irg 



N sin (w(m - n + g) /N) 



(6) 



where $ m are independent random phases uniformly dis- 
tributed between and 2tt and g is some fixed constant. 
A closely related model, introduced in [l8| and investi- 
gated in corresponds to a family of quantum maps 
obtained by quantization of a classical map on the torus. 
These quantum maps have exactly the form (|6]) , but with 
g now depending on the matrix size, namely g = aN with 
a fixed. Spectral statistics of these ensembles have been 
shown to be of intermediate type [2(| • 

We will also consider the well-studied CrBRM ensem- 
ble, which is the set of N x N symmetric or Hermitian ma- 
trices whose elements M mn are independently distributed 
Gaussian variables with zero mean and variance 



\M nn \ 2 ) 



\M„ 



1 



9 



m 7^ n. (7) 



Here /3 = 1 if the matrices considered are real symmetric, 
or /3 = 2 if they are complex Hermitian. 
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III. MULTIFRACTAL DIMENSIONS AND 
SYMMETRY RELATION 

We want to investigate multifractality properties of 
the above ensembles of random matrices. Multifractality 
manifests itself in a non-trivial scaling of eigenfunction 
moments. Let *f/ n (a) be a normalized eigenfunction of 
the matrix M mn corresponding to eigenvalue X a , so that 

N 

M mn V n {a) = \ a V m {a). (8) 

71=1 

Mean eigenfunction moments are defined for any real q 
as I q = (J^j \ 1 &j\ 2q ), where the average is taken over all 
random variables of the model and over a small window 
of eigenvalues. More precisely, we set 

1 n 1 N 

a=l 1 K ' j=l 

where p(E) is the total mean eigenvalue density 
1 N 

p(e) = ^(J2 s ( e - x ^)- ( 10 ) 

a=l 

Multifractal dimensions D q are defined as the scaling ex- 
ponents in the large N behavior of these moments, 

I q ~ Ar-A,(<2-i) _ (11) 

N— >oo 

In our matrix models, the dimension d of the underlying 
space is 1, so that anomalous dimensions are given by 
A q = (D q — l)(q—l). Both analytical and numerical cal- 
culation of these exponents is a difficult task. Indeed, it 
is so far not possible to have direct access to exact analyt- 
ical expressions, and numerical calculations require diag- 
onalization of large matrices. There is however a way to 
get some insight into analytical expressions by construct- 
ing perturbation series for the eigenfunction moments for 
small or large values of the coupling constant [5j, |g, l21| . 
For moments of order q with q > 1/2, the lowest order 
of perturbation series at small and large g has a univer- 
sal dependence on q for all critical ensembles considered 
above. In the regime of weak multifractality, the unper- 
turbed eigenvectors are delocalized, and at first order 

D q ~ 1 - tq, (12) 

where the constant t depends on the system. It is equal 
to l/(27r/90) for CrBRME, and for CM models, where 
this regime corresponds to g — > oo. For RS ensemble, 
weak multifractality regime is reached whenever g is close 
to a nonzero integer, and t = (g — k) 2 /k 2 , where k is 
the integer closest to g [l2| . Equation (TP2")) is valid for 
t\q\ <C 1. In the regime of strong multifractality (|<7| <C 
1), the unperturbed matrices (corresponding to g = 0) 



are diagonal and their eigenvectors are localized. The 
first-order correction is given by 

Dq g ^\g\ P {E)s ^ , (13) 

again valid only for q > 1/2. Here p(E) is the density 
of states of unperturbed system, which, at that order, is 
equal to the actual density of states (jTUJ) . The constant 
s in (fT3"|) can be obtained as the coefficient in front of the 
logarithmic term 

1 N 

(]vEiM^ 2skJV ' (14) 

The explicit expressions of s can be found e.g. in [l2j]. 
It is equal to 1 for RS, CM r and CM^ ensembles, I/a/tt 
for symmetric CrBRME, i/ttJS, for Hermitian CrBRME, 
and [h/tv] for CM t , where [.] denotes the integer part. 

In the weak multifractality case at leading order, the 
fact that D q is linear in q (Eq. (TT2"|) ) implies that the 
symmetry relation holds. Indeed, A g ~ — q(q — l)t is au- 
tomatically symmetric with respect to the substitution 
q 1 — q. In [22[ it was demonstrated analytically that 
for CrBRME and ultrametric ensembles the second or- 
der term also obeys this symmetry. In the same paper 
the third order term has been calculated numerically for 
these models and it has been found that it also fulfilled 
Results in [16j suggest that these results extend to all 
orders of perturbation theory in the weak multifractality 
regime. 

The situation in the strong multifractality limit is less 
clear, since (fTB")) is valid only when q > 1/2 5]. The 
purpose of the present paper is to obtain an expression 
for Dq in this regime when q < 1/2. We show that it 
can be derived within the standard perturbation series 
approach, and yields, for q < 1/2, 



2a — 1 v " 1 1 2 v 

D q ^ n ^ + M9\P(E) s - — -1 . (15) 
5^0 q - 1 ( q - l)T(-q) 

The existence of a g-independent term in (fT5"|) and its 
absence in (TTB")) when q > 1/2 looks, at first glance, as 
the necessity of non-analytical contributions. In fact, as 
we shall see, the perturbation series is obtained in an 
even simpler way than the one used to derive (|13[) . 

A remarkable feature of (fT5|) is that it coincides with 
what would be obtained from (|13l) assuming that the 
symmetry ([3]) holds. Therefore this perturbation theory 
approach shows that the symmetry ([3]) is correct at order 
zero and at order 1 in the coupling constant g for all 
ensembles considered. 
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IV. ZERO-ORDER TERM 

Let us consider a set of random Hermitian matrices of 
the form 



: Pm $mn 9\\ ^ ' ] mn) •mn j 



(16) 



where p m are independent identically distributed random 
variables distributed according to some probability dis- 
tribution a{p). The off-diagonal elements [i mn may con- 
tain random variables, but all are independent on the p m . 
Matrices for CM and CrBRM ensemble are obviously of 
that form. The case of the RS ensemble will be discussed 
separately in section IV B[ 

Standard perturbation theory [23j applied to matrix 
(dn]) for g <C 1 yields eigenvalues 

\ a = Pa + g 2 Y,^- + °(9 3 ) (17) 



Pa -Pn 



and eigenfunctions 



E 



l^mnf-^na 



Pn ~ P, n ^ (Pa - Pn)(Pa ~ Pm) + ^^ 

n^m,a. 

(18) 

for m 7^ a, and ^ a (a) = 1 + 0(g 2 ). In the first non- 
vanishing order, for fixed a, we neglect the term in g 2 in 
(0 and get 

AT 

J2\^M\ 2q S(E~X a ))= (19) 



m^a 



since by assumption the p 3 are independent on the /x mn 
and A Q = p a + 0(g 2 ). In order to calculate the averages 
one has to integrate over random variables p n distributed 
according to the law a(p). The average over variables p a 
and pm with m ^ a reads 

S(E- Pa ) \ f S(E- Pa ) 
= / ap a dp m a{p a )a(Pm) 



(Pa - Pm) 2q 



= o-(E)A, 



(Pa - Pm) 2q 

(20) 



where we have set 



a(p) 

A « = I dp jE—p^- 



(21) 



This integral is finite whenever g < 1/2. From Eqs. (flU)) 
and (fT7l) one has 



p(E) = a(E)+(D(g 2 ), 
thus the moment © reads 



^ = 1 



N 



9 EE<i^i 29 >- 



(22) 



(23) 



The summation in (|23p can be performed in the same 
way as was done e.g. in [l2f to obtain fH|l . but instead 
of the IniV behavior of (Tl4| . for q < 1/2 these quantities 
will scale as TV 1-29 when AT — > 00. For instance in the 
case of CM r model, [i mn = i/(m — n), and for q < 1/2 



N-l 



1 1 _ 2 N — k 



N £l, \ m ~ n \ 2q N _ k2q AT-Kx, (l-2gr)(l-qr)' 



m,n 



so that for large AT, the mean moment J g scales as 



N- 



00 9 (l-2q)(l-q) 



2 <i. 



TV 1 - 2 ?. 



In the trigonometric case, one has 



1 

N 



E 

m,n 



2-/ 



N 2 i t 



,2g 



(/i(m — n)/N) 



(24) 



(25) 



(26) 



- 2 M 2<? iV 1 ~ 2 ' ? / dy- , 

io sin 2 «( ra )' 



1 



?7 



so that 



I ~ ffV%/ dy^-% 1 ^. (27) 

For the other ensembles considered, the calculation is 
similar (technical details are the same as the ones encoun- 
tered for calculating (|T4")) . see [Hj]). All these ensembles 
yield the same N 1 ~ 2q behavior. 

The definition (fTT|) of fractal dimensions then allows 
us to extract the leading order term for multifractal di- 
mensions in the domain q < 1/2 as 

2q- 1 



D„ 



<1 



1 



(28) 



Note that for q > 1/2, the various integrals appearing 
in the above calculation do not converge. Thus, this 
usual perturbation theory approach breaks down. Dif- 
ferent techniques were used in the case q > 1/2, yielding 
an expression for D q given by Eq. (|1 3|) . At order zero 
this expression (|13l) gives D q = 0. In particular, together 
with Eq. (|28p. these results show that the symmetry re- 
lation ([3]) is fulfilled at order zero of perturbation theory. 

In order to assess the correctness of the g 2q behav- 
ior in the leading-order term in the expressions (|25[) 
(j27l) , we take as an example the case of the trigonometric 
Calogero-Moser model in Fig. [TJ According to (|27p. the 
ratio I q /N 1 ~ 2q is expected to behave as C q g 2q for small 
values of <?, with C q a factor depending on q. In log- 
arithmic scale we thus expect straight lines with slopes 
2q. The comparison with this predicted power, given in 
Fig. [TJ is very good. 



V. FIRST-ORDER TERM 

To get the first-order correction to the above result 
one has to take into account the term in g 2 in (Tj"5|) . The 
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FIG. 1. (Color online) Circles: Slopes a obtained by a linear 
fit In^q/N 1 - 2 ") = a\ng + b, versus q, for CMt ensemble. The 
slope is extracted from values of g from 5.10 -5 to 2.10 -3 . 
Matrix size is N = 2 10 + 1, average is taken over 64 central 
eigenvectors (i.e. around E = 0) and 1600 realizations. Solid 
line is the best linear fit a — —0.186 + 1.99q (to be compared 
with the expected a = 2q+constant from Q27| l1, 



modulus square of the wave function with m ^ a is then 
given by 



In this case (|29|) gives 
\^m(a)\ 2q = 

9 



9 2q \^ ma \ 2q 

(Pa -Pm) 2q 
2 



(30) 



x 1 + 



■ E 

n,n' ^m,a 



(j>a -Pn)(Pa ~ Pn') 



Note that at this order it is still possible to replace a(E) 
by p(E) and 5{E — X a ) by 8{E — p a ) since corrections 
are of higher order. The average (S(E — p a )/ (p a — p m ) 2q ) 
over variables p n in (|19p is now replaced by 



1- 



E 



S{E- Pa ) 

(Pa - Prn) 2q ^ ' \^ma\ 2 (Pa _ Pn)(Pa - Pn') 



(31) 

Since the variables p n and p n i appearing in (|31[) are dif- 
ferent from the variables p a and p m , the average (f3"Tj) can 
be done separately over variables p n ,p n i. For the term 
involving these variables, the correction to the order zero 
is given by 



9 



E 



-l) , (32) 



where the average is performed only over variables p n 
with n ^ m,a. For fixed the main contribution to 
(|32p . upon integration over random variables, will come 
from terms with p n ~ p Q . The leading-order correction 
(|32| can in fact be replaced by 



l* m («)| 5 



C.C. 



(Pa - Pm) 2 

f J -mri^nafJ' rna 
Pa-Pn 



^2 MmnMnaM mn 'M n ' a 
77* (Pa -Pn)(Pa ~ Pn') 



(29) 



(c.c. denotes the complex conjugate). As for all critical 
systems, the dominant correction in such formulae when 
g — > corresponds to transitions between two states n 
and a with a small energy difference p n — p a of the order 
°f 9 0- Taking more states in resonance into account 
gives higher-order corrections @ . Calculation of the first- 
order term depends on whether the off-diagonal matrix 
elements are pure imaginary or not, since in the pure 
imaginary case the term in g 3 vanishes in (|29l) . Since 
matrices for CM ensembles have this feature, we first 
consider this case. 



A. CM ensembles 

CM ensembles defined by Q have the property that 
their off-diagonal elements are pure imaginary numbers. 



E 



l 



9 



l^mal 2 (Pa ~ Pn)' 



1 



(33) 



In order to show this, one has to check that the difference 
between (|32|) and (f33|) is of higher order in g. This can be 
done by upper bounding the integrand of the difference 
(in the integral over p n ) uniformly in g by a function 
whose integral converges at g = 0. Note that a similar 
problem occurs in the calculation of the zero-order term 
(|20[) which we obtained in the previous section: for fixed 
g the average (f31"j) diverges at p n — p a , but nevertheless 
(f3"2"| converges to when g — > 0, which proves that (|2"0|) 
gives the correct zero-order contribution. 

The term (|33[) is a sum over n. For a given fixed index 
n, the contribution is 



Jn(9) = / <lPnCr(Pn) 



1 



9 



| P'mn^na \ 



\Vma\ 2 {Pa - Pn) 2 

Changing variables by the transformation 



1 



Pn 



Pa 



\f^mnf-^na . 

9—, r- * 

\f^ma | 



(34) 



(35) 



and keeping only lowest-order terms in g in the expansion 
of (j{p n ), the contribution (p4|) gives 



Jn(g) 



\9\—, ^cr(Pa) 



i v 

1+ ^> 1 



dt . 

(36) 
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The integral can be calculated for q < 1/2 (see e.g. [2J]), 
yielding 



B q = 



1 \9 



d* = -- 



^(1/2 -g) 

r(-g) 



(37) 

The sum of all contributions in (|33p gives the first-order 
correction term 



S(E-p a ) 

(Pa - Prn) 2q 



E 



| f^mnf^na 



a(p a )B q 



(38) 



I l^mnl^noi \ 
\f^"ma | 



using the result of (|20p. Putting together equations ©, 
(fT9l) , (l3"T1) and (l38l) , the first-order correction to moments 



I q is 



A^ 



|2q Iftmnftna \ 
\l^moi\ 



(39) 

In the case of CM r model where /z mn = i/(m — n), we 
need the asymptotic large- N behavior of 



-y. 



1 



(to — a) 



2g 



E 



l 



l 



to — n a — n 



(40) 



This sum is dominated by two regions n = to + k with 
1 -C fc *C to and n = a + fc with 1 <C k <C a. In both 
regions the sum over k diverges logarithmically, so that 



E 



l 



l 



to — n a — n 



E 



i 



21niV 



| TO — 71 1 



E 



l 



\a — ra | 



(41) 



Using (j2~4")l and (|4*T1) for these two regions, we conclude 
that the behavior of the sum (gDJ is - 47V 1_29 lniV 
for large AT. As we show in the following sections, for 
other ensembles similar calculations yield a behavior in 
4sN 1 ~ 2q \nN, where s is the same as in (fT4")l . Thus at 
first order 



(42) 



where B q is given by (f3"T|) . This can be seen as a first- 
order expansion of C q (g)N^ D ^Kq~i) at small values of 
g. We thus get the first-order correction to the fractal 
dimensions as 



Vl + t | )W «, ,43, 

9-1 (? - i)r(-<?) 



which coincides with Eq. (|15l) . As already mentioned, 
together with the earlier result ([13"]) . this confirms the 
symmetry relation Q at first order of perturbation series 
at small g for models 



B. RS ensemble 



For RS ensemble, expanding ([5]) at small g yields, up 
to second order in g, 



M mm ~ e 1 *™ 1 - 



7rV(iV 2 -l) 
6N 2 



(44) 



for diagonal elements and 



(45) 



sin (7t(to — n)/N) 

ng cos (7r(m — n) /AT) \ 
~ A^sin 2 {■K{m-n)/N)j 



for off-diagonal elements with to ^ n. Thus at that order 
the matrix M takes the form (|16[) if we set p m = M mm 
and M m „ = Keeping only terms up to order # , 

the expression (fig]) now gives 



*m(a) = 



Try 



AT e 1 ' 



e 1 " 



1 



sin(7r(TO — a)/N) 



(46) 



ixg cos (7t(to — a)/N) -ng 
N sin 2 (7t(to — a)/N) + iV 



^ y 

AT P i4 



sin (7t(to — n) /N) sin (7r(n — a) /N) 
Using the identity for m a 
1 



iV 



V - 



(47) 



^ sin(7r(TO — n)/N) sin(7r(n — a)/N) 

2 cos (n(m — a)/N) 
sin 2 (7t(to — a) /AT) 

and the fact that 

e'*™ i eft _ <J> i 

= -cot " a -~, (48) 
2 2 2' V ' 



one can check that two terms cancel in (|46l) . so that 
^ m (a) reduces to 



*m(a) 



^ sir 



7r.g e 



AT e 1 *" - c 1 ' 



sin(7r(m — a)/N) 



(49) 



icot((*„-$„)/2) 



2N ^ sia Mm - n) /N) sin Mn- a)/ N) 



This has a form similar to JT8J), with a pure imaginary 
first-order correction term, for leading-order terms where 
&a ~ ^n- K can be treated essentially in the same way, 
yielding the same expression (|43[) for D q , with s = 1. 
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C. Critical random matrix ensemble 

In the above derivation we used the fact that for models 
(gD and © the term in g 3 in (gSJ vanishes. For CrBRME 
(O this is not true and this case requires a different treat- 
ment. One can consider this ensemble as a particular case 
of (TlBl) with the substitution = M^k and (at first order 
in g) n mn distributed as independent Gaussian variables 
with zero mean and variance 



(| Mr, 



1 



(m — n) 2 



(50) 



Since pL ma and (fTHj) are independent complex 

Gaussian random variables, ^ m (a) appears as a sum of 
random variables whose variances are given by (|50p . It 
is thus itself a complex Gaussian random variable Z of 
mean and variance 

{\Z\ 2 ) = 9 \ , 2 ( (|Mma| 2 )+5 2 E <lMmn|2>lMnQl 



(Pa -Pra) 2 



(Pa -Pfe) 2 



(51) 

(at this stage we consider the fi na fixed). For q > — 1, 
the mean value of the gth power of Z can be calculated 
as 



(\z\*) = (\z\ 2 )*r( q + i) . 



so that 



(\* m (a)\^)=T( q+ l) g q ^ mal I" 

\Pa Prn ) 



(52) 
(53) 



1 + E 



(|Mrrcn| )lMria| 
(|Mm Q | 2 )(Pa - Vnf 



The treatment is thus the same as in the previous section 
(compare with ([3D])), apart from an overall factor T(q+1) 
and with the replacement 



\2q 



(I Mr. 




(54) 

In view of (|5T)|) the calculation is exactly the same as for 
the CM r ensembles and yields (|4"3")l for _D g with s = 1. For 
q < — 1 the mean value of the gth power of Z is infinite, 
and mean quantities should be normalized by the typical 
value of Z, as was discussed in detail in [25J. Again, 
the symmetry ([3]) is verified at this order of perturbation 
theory. 



VI. NUMERICAL RESULTS 

In previous sections we showed that the symmetry re- 
lation d3j holds at lowest orders of perturbation the- 
ory for small coupling constant g. For large values of 
the coupling constant (weak multifractality limit), re- 
sults have been presented in pjj]. In order to obtain 



results for general coupling constant (between the two ex- 
treme regimes where perturbation theory holds), we have 
to resort to numerical investigations. Here we present 
numerical results for the matrix ensembles considered 
above at various values of the coupling constant. Ma- 
trices are diagonalized and moments of the eigenvectors 
are calculated for different matrix sizes. The large-iV 
behavior of moments is then extracted via a fit of the 
form log(^ i |\&i| 2q ) = a + b\ogN + c/N. For negative 
values of q, moments are replaced by a local averaging 
l*4i+i| 2 ) 9 to avoid divergences. 

Figures [5] and 0] display the multifractal dimensions 
D q as a function of q for CM r and CM t ensembles re- 
spectively. As CM matrices are Hermitian with a den- 
sity p(E) depending on the energy window considered, 
the averages are performed only over eigenvectors whose 
eigenvalues are in the vicinity of E = (i.e. in the center 
of the spectrum). For RS ensemble, multifractal dimen- 
sions D q are shown in Fig. [5] In this case, matrices are 
unitary, and averages are performed over all eigenvectors. 
Figure [5] displays multifractal dimensions of CrBRME. 
For all these models, the plots of D q follow Eqs. (fT3"|) 
and (fTS")) for small values of the coupling constant. This 
was in fact checked indirectly in where plots of D q 
for small g and q > 1/2 were shown to agree with (fT3|) . 
while plots for q < 1/2 agreed with the expression ob- 
tained from (|13[) assuming the symmetry ([3J). Since this 
symmetry is now proven analytically in the small g limit, 
the plots for q < 1/2 indeed coincide with (|T5|) . 

Figures El El El and M are a direct check of the sym- 
metry of anomalous dimensions A q for ensembles CM r , 
CM t , RS, and CrBRM, respectively. As expected from 
the results obtained above, the symmetry is well verified 
in both limits of small and large coupling constant (or 
in the case where a is close to an integer for RS ensem- 
ble). Our figures show that the symmetry is also quite 
accurately verified for all other values of the coupling 
constant, both for q £ [0,1] (see insets of figures HHPJ , 
and for larger values of q. 



VII. INTERMEDIATE MAP 

Although the symmetry of anomalous exponents can 
be observed in a variety of systems, there are however 
situations where no symmetry can be observed in the case 
of strong multifractality. This is the case for instance for 
the ensemble of N x N unitary matrices given by 



M - e 1 *" 



sin(7raA) 



N sin (7r(m - n + aN) /N) ' 



(55) 



where $ m are independent uniformly distributed random 
variables, and a some fixed number. This ensemble is 
identical to RS ensemble, but with a constant g = aN 
which now depends on the size of the matrix. Spectral 
statistics of this model turn out to be very different from 
those of the RS ensemble [2(|. For rational a — m/b it 
displays intermediate spectral properties which depend 
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FIG. 2. (Color online) Fractal dimensions D q as a function 
of q for CM r ensemble for (from bottom to top at q = 4) 
3 = 0.005 (red), 0.025 (green), 0.05 (blue), 0.15 (violet), 
0.25 (brown), 0.4 (orange). The pk are independent random 
variables distributed according to a Gaussian with mean 
and variance 1. Matrix sizes for numerical fit are TV — 2 n , 
8 < n < 12. Average is performed over the TV/16 eigenvectors 
closest to the eigenvalue E — 0. Number of random realiza- 
tions of the matrix is between 2560 for N — 2 s and 80 for 
2 12 ; error bars correspond to the standard deviation. Dashed 
curve is the function q \-¥ (2q — l)/(g — 1). 



T 




I I I I i L 

-3-2 2 

q 



FIG. 4. (Color online) Fractal dimensions D q as a function 
of q for CM t ensemble for /i = 2n/N and (from bottom to 
top at q = 4) g = 0.005 (red), 0.025 (green), 0.05 (blue), 
0.1 (violet), 0.15 (orange), 0.25 (indigo), 0.4 (brown). The 
Pk are independent random variables distributed according 
to a Gaussian with mean and variance 1. Matrix sizes for 
numerical fit are TV = 2 n , 8 < n < 12. Average is performed 
over the TV/16 eigenvectors closest to the eigenvalue E = 0. 
Number of random realizations of the matrix is between 2560 
for TV = 2 8 and 80 for 2 12 . Dashed curve is the function 
9 rt(2?-l)/( 9 -l). 
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FIG. 3. (Color online) Anomalous dimensions A q — (D q — 
l)(q— 1) (solid) and Ai- q (dashed) as a function of q for CM r 
ensemble, same data and same color code as in Fig. [2] g in- 
creases from bottom curve to top curve at q = 4. Inset: same 
data zoomed in, circles correspond to A q , stars to Ai_ 9 . For 
almost all data points stars lie inside circles, which indicates 
symmetry of A q . 



FIG. 5. (Color online) Anomalous dimensions A q (solid) and 
Ai-q (dashed) as a function of q for CMt ensemble, same data 
and same color code as in Fig. [4] g increases from bottom 
curve to top curve at q = 4. Inset: same data zoomed in, 
circles correspond to A q , stars to Ai_ 9 . 



on the remainder of mN modulo b. The mult if r act ality 
spectrum for such an ensemble was investigated in |17| . 
The strong multifractality regime is illustrated in Fig.fTUl 
for rational values a = 1/b with b = 3, 5 and 9. As 
can be seen, its features are quite different from those 
of other ensembles: for the smallest denominators the 



multifractal dimensions are quite far from the limiting 
value (2q — l)/(q — 1). Note that in this model there is 
no clear way to extract a perturbation series expansion. 
As for the symmetry of anomalous dimensions, it does 
not hold in general. Indeed, one can see from Fig. [11] 
that the value of Ai_ q clearly deviates from A q for small 
values of the denominator b. 
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FIG. 6. (Color online) Fractal dimensions D q as a function of 
q for RS ensemble for (from bottom to top at q — 4) a = 0.01 
(black), 0.05 (red), 0.1 (green), 0.2 (blue), 0.3 (violet), 0.5 
(brown), 0.7 (maroon), 0.9 (cyan). Matrix sizes for numerical 
fit are N = 2 n , 8 < n < 12. Average is performed over all 
eigenvectors. Number of random realizations of the matrix is 
between 1024 for N = 2 s and 64 for 2 12 . Dashed curve is the 
function q H> (2g — !)/(</ — 1). 
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FIG. 7. (Color online) Anomalous dimensions A q (solid) and 
Ai_, (dashed) as a function of q for RS ensemble, same data 
and same color code as in Fig. [6] g increases from bottom 
curve to top curve on the right; for better clarity, data cor- 
responding to a = 0.05 and 0.1 have been removed. For RS 
ensemble, the plot for Ai_ 9 is almost indistinguishable from 
the plot for A 9 . Inset: same data zoomed in, circles corre- 



FIG. 8. (Color online) Fractal dimensions D q as a function of 
q for CrBRME with /3 = 2. Constant is (from top to bottom 
at q = 4) 1/g = 0.1 (black), 0.5 (red), 2 (green), 10 (blue), 30 
(violet), 60 (orange), 80 (indigo), 100 (brown). Matrix sizes 
for numerical fit are N — 2 n , 8 < n < 12. Average is per- 
formed over the N/16 eigenvectors closest to the eigenvalue 
E = 0. Number of random realizations of the matrix is be- 
tween 2560 for N = 2 8 and 80 for 2 12 . Dashed curve is the 
function q (2q — l)/(q — 1). 




spond to A q , stars to Ai_ 



FIG. 9. (Color online) Anomalous dimensions A q (solid) and 
Ai_ 9 (dashed) as a function of q for CrBRME, same data as 
in Fig. [8] (g increases from bottom curve to top curve on 
the right; for better clarity, data corresponding to 1/g = 30, 
60 and 80 have been removed). Inset: same data zoomed in, 
circles correspond to A q , stars to Ai_ 9 . 



VIII. AN EXACTLY SOLVABLE CASE 

We now briefly turn to an example closely related to 
the above RS ensemble, and where multifractal dimen- 
sions can be calculated exactly. Let us consider vectors 
of the form 

_ sin(Tra) 
p ~ Nsm(7r(p + a)/N) [ ' 



for 1 < p < N, a being a fixed real number. These 
vectors correspond to column vectors of the matrix ([5]) 
for RS ensemble. Namely, M^ p = c'^^-p, and since 
Mfcp is unitary, vectors (|56[) are normalized by 

N-i 

p=0 
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where 



FIG. 10. (Color online) Fractal dimensions D q as a function 
of q for intermediate map for a = 1/fe with (from bottom to 
top on the right) b — 3 (black), 5 (red), 9 (blue). Matrix 
sizes for numerical fit are the smallest integer N' > N — 2", 
8 < n < 12 such that N' = 1 [6]. Average is performed over 
all eigenvectors. Number of random realizations of the matrix 
is between 1024 for N = 2 s and 64 for 2 12 . Dashed curve is 
the function q H> (2q — l)/(q — 1). 
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FIG. 11. (Color online) Anomalous dimensions A q (solid) 
and Ai_ 9 (dashed) as a function of q for the intermediate 
ensemble (|55[) . same data as in Fig. 1101 Inset: same data 
zoomed in. 



Our goal in this section is thus to calculate fractal dimen- 
sions of a system whose eigenvectors would be columns of 
the unitary matrix Mk P - As we now show, multifractal 
dimensions for the vectors (|56l) coincide with the zero- 
order term in the perturbation series expansion of mul- 
tifractal dimensions for eigenvectors of CM, RS and Cr- 
BRM ensembles in the strong multifractality regime, that 
is, Eqs. flT3J) and $T5§ when g = 0. 

Moments of the vectors ([56]) are given by 



N-l 
p=0 



|T l2o /sin(7ra)\ 2 9 



(58) 



N-l 

f N (q,a) = Y 

p=0 



i{ir(p + a)/N)\ 2q 



(59) 



In order to get fractal dimensions, we need to calculate 
the asymptotic behavior of this sum in the limit of large 
N. The usual way of estimating such sums is the Eulcr- 
Maclaurin formula. Here we need only the leading term 
in powers of N. 

The sum (|59p is periodic in a with period 1 and it is 
sufficient to consider < a < 1. Let M = [N/2] and 



M-l 



1 



^ sm^(7r(p + a)/N) 



(60) 



In the case where q < 1/2, the sum gtq{q,a)/N can be 
easily lower and upper bounded by Riemann sums which 
for large N both converge to the integral 



1 



7T/2 



dt 



r(i/2- g ) 

sin 2 »(i) 2VH\1 - q) 



(61) 



For q > 1/2, using the inequality x — x 3 /6 < sins < x 
over [0,7r/2] one gets 



M-l 

E 

p=0 



1 



< 



7T q 9N(q,a) 



(p + a) 2q 

ir 2q g N (q,a) 
N 2 i 



N 2q 

M-l 



(62) 



1 



1 



" ^ {v + a) 2q [l-^{p + a) 2 ] 2q 



p= 



The left-hand side sum tend to £(2g, a) when N — > oo, 
where £ is the Hurwitz zeta function (see e.g. [Hj]) denned 
when s > 1 as 



1 



p=0 



(p + a) s 



(63) 



It is not difficult to show that the right-hand side also 
tends to the same limit when q > 1/2. Thus from (|6"2"|) 
we obtain that ir 2q gN(q,a)/N 2q converges to ((2q,a). 
Summing up the results we get for q > 1/2 

J im -^fN(q,a) = ^-k(2q,a) + a2q,l-a)] . (64) 
and for q < 1/2 

1 r(l/2-g) 
lim —fN(q,a) = . 65) 

Af->oo A y'TT 1 (1 — q) 

For q = 1/2 one has an additional log N term 



lim — /jvd,a) — ► § 



lnA r -ln^--(*(a)+*(l-a)) 



where ^/(z) is the digamma function. 



(66) 
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(small coupling constant). However in the latter case, 
only the region q > 1/2 had been investigated, since the 
methods used for this region cannot be directly applied 
to the region q < 1/2. 

In this paper we filled this gap by developing a per- 
turbation theory valid in the region q < 1/2 for small 
coupling constant. We considered critical random ma- 
trix ensembles built from Lax pairs of classical integrable 
iV-particle systems, as well as the well-studied critical 
banded random matrix ensemble. We derived a pertur- 
bation expansion of the moments of eigenfunctions 



|^| 2 9) =C q (g)N- D ^ q -^ 



(68) 



FIG. 12. (Color online) Fractal dimensions D q as a function 
of q for column vectors of Ruijsenaars matrix (vectors (|56[0 
with a — 0.5, fit over sizes 2 s to 2 12 . Dashed: theoretical 



expression D q 
q > 1/2. 



(2q - l)/(q - 1) for q < 1/2, D q = for 



and obtained the leading terms of the perturbation series 
in g for multifractal dimensions, 



D q (g) 



2q-l 



A\g\p{E)s 



(?-i)r(-«) 



(69) 



From these results it follows that the fractal dimensions 
T)6"|) is given by 



D q of vectors 







for q > 1 /2 



D «=< *m for g <l/2 

? - 1 



(67) 



which exactly obeys the symmetry relation ([3]) and coin- 
cides with the perturbation series result at g = (dashed 
line in Figs. [U0] and E]). ^ n Fig- [HI we present the plot of 
D q for vectors ([56]) with a — .5. Numerical results agree 
very well with our analytical result (|67p when g > 1 and 
q < 0. However when q gets close to 1/2, both terms 
in and (l6"5j) contribute for finite values of N and 
it is difficult to numerically recover Eq. (|67p due to the 
existence of the logarithmic term in (|66p . 



IX. CONCLUSION 

We considered critical random matrix ensembles whose 
elements decrease as inverse powers of their distance to 
the diagonal. It is known that for such ensembles eigen- 
functions have multifractal properties, which are charac- 
terized by a nontrivial set of multifractal dimensions D q . 
Perturbation theory, which is the main analytical tool to 
access these multifractal dimensions, was successfully ap- 
plied before, both in the weak multifractality limit (large 
coupling constant) and in the strong multifractality limit 



At order zero, C q ~ g 2q and multifractal dimensions are 
given by a simple but nontrivial term D q = (2q— l)/(g — 
1). The next order is obtained by considering the second- 
order term of perturbation theory. Higher-order terms 
can be obtained straightforwardly by the same method. 

It turns out that the above formula for D q at q < 1/2 
exactly coincides with the expression obtained from D q 
for q > 1/2 when assuming that a symmetry holds 



for anomalous dimensions A 



Ai_ g , where A q 



(D q — l)(q — 1). Therefore, our results also prove that 
this conjectured symmetry relation is valid for the mod- 
els considered at leading orders of perturbation theory 
in the strong multifractality regime. We also performed 
detailed numerical investigations which indicate that this 
symmetry property is in fact fulfilled for all coupling con- 
stants. 

Nevertheless, this symmetry is not an intrinsic prop- 
erty of all critical random matrix ensembles. Indeed, by 
a small modification of one of our models, we constructed 
an ensemble where numerics clearly show that such a 
symmetry is absent. Therefore the range of validity of 
this symmetry relation remains an open problem (see 
however the recent paper (lrjj). 
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